Extrapolation method in shell model calculations with deformed basis 
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An extrapolation method in shell model calculations with deformed basis is presented, which uses 
a scaling property of energy and energy variance for a series of systematically approximated wave 
functions to the true one. Such approximated wave functions are given by variation-after-projection 
method concerning the full angular momentum projection. This extrapolation needs energy variance, 
which amounts to the calculation of expectation value of square of Hamiltonian H^. We present the 
method to evaluate this matrix element and show that large reduction of its numerical computation 
can be done by taking an advantage of time-reversal symmetry. The numerical tests are presented 
for fp shell calculations with a realistic residual interaction. 

PACS numbers: 21.60.Cs 



Extrapolation is an intriguing and useful concept in 
computational physics in quantum many-body problems 
in condensed matter physics, nuclear structure physics 
and other fields. The Hilbert spaces of quantum many- 
body systems can be finite but essentially large or infi- 
nite. To handle such Hilbert spaces completely by nu- 
merical methods, there exist various difficulties. How- 
ever if we can handle its dominant but small subspace, 
the contribution of the rest huge subspace can be taken 
into account by extrapolation. More concretely we can 
estimate the exact energy by extrapolation if we know a 
scaling property how to change energy as a function of a 
certain physical quantity. 

In the nuclear structure physics, shell model is one of 
fundamental frameworks. For shell model calculations, 
one needs handle very huge Hilbert space. One of conven- 
tional approaches is diagonalization with spherical single 
particle basis states, which recently have been applicable 
to quite large-scale problem with dimension up to 10^0. 
However, for further larger spaces, diagonalization is im- 
possible. To overcome such a limitation of diagonaliza- 
tion, various methods 

have been proposed and have been developed. Among 
them, extrapolation methods 0, H, in III have 
been developed, combined with diagonalization of the 
Hamiltonian matrix evaluated with spherical basis states. 
For instance, exponential convergence method (ECM) % 
can predict the exact energy by assuming exponential 
behavior of energy as a function of the number of ba- 
sis states and empirical successes of its assumption have 
been reported with various truncation schemes jSj, iQ, ilO] . 

We have also proposed an extrapolation method in 
shell model calculations, combined with Lanczos diago- 
nalization for a series of truncated spaces 0, , where 
we take an advantage of a scaling property between the 
energy difference SE and energy variance AE for a series 
of systematically approximated wave functions. Such a 
scaling property has been originally introduced [T^ ITsI 
in the study of Hubbard models. We have introduced 
this scaling property into shell model calculations [l^ . 
Moreover, we have shown that, for precise estimation, 
higher order term of expansion can be taken into consid- 



eration .13 ] . Very recently this scaling property has been 
successfully applied to the no core shell model [lTI| . 

Thus two kinds of extrapolation methods d, 0, 0, 
[m Il2l IT^ have been proposed and are quite promis- 
ing. However they have a common drawback that they 
strongly depend on the diagonalization with spherical ba- 
sis states. The dimension of necessary truncated sub- 
space heavily depends on each nucleus as the dimension 
of whole shell model space does. 

In the present paper, we propose an extrapolation 
method in the shell model calculations with deformed 
bases, which are expected to be suited for deformation 
of nuclei. As our extrapolation scheme depends on only 
energy and energy variance, we can apply it to deformed 
basis. We will show that the exact energy can be ex- 
trapolated by the behavior of energy evaluated by single 
Slater determinant with angular momentum projection. 

First we briefly summarize our extrapolation scheme. 
We consider a series of systematically approximated wave 
functions as \'ip2),'''- Here we assume that they 

are systematically generated and their energies becomes 
lower as El > E2 > ■ ■ ■ where Ei is an energy of 
Later we will discuss how to prepare them in deformed 
basis in detail. The approximated wave function contains 
large fraction of true wave function and contamination 
of rest wave function is small. By expanding this small 
value, a following relation is obtained. We define the 
difference 6E between the lowest energy eigenvalue (H) of 
approximated wave functions and true energy eigenvalue 
{H)o as SE = (H) — {H)o- The energy variance AE 
of approximated wave function is also defined as AE = 



As a first order approximation, we can show 



a proportionality as 5E oc AE by expanding 5E 
as a function of AE. We can take into account second 
order effects for precise estimation . By these scaling 
properties, we can estimate the exact energy by taking 
AE 0. 

Next we will discuss how to prepare such systemat- 
ically approximated wave functions in deformed basis. 
We consider an angular momentum projected deformed 
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wave function such as 

\i^.LM)=Y.gKPLK\^) (1) 
K 

where P^jj^ is angular momentum projection operator 
and g's are coefficients. The wave function \^) is a Slater 
determinant as = 11 "^a 1^) where |0) is a vacuum state 
and ajj's are creation operators on deformed orbit a and 
are defined as a^. = "^DaAcj. The c| denotes creation 

i 

operator of spherical orbit i and D's are coefficients and 
are determined by the following variational equation: 

[ {WJ,M IVJM) J 

In order to solve above equation, we minimize the en- 
ergy expectation value with angular momentum pro- 
jected wave function concerning Da^i [l^ . This can be 
achieved by gradient method where we evaluate the gra- 
dient vector in the projected energy surface and then we 
change the wave function along the steepest descent line. 
By iterating this procedure, the variation of the angu- 
lar momentum projected energy is achieved. In general, 
obtained wave function in this procedure corresponds to 
local minimum in angular momentum projected energy 
surface. By examining several local minima, we can ob- 
tain lowest minimum in the practical calculations. Here 
we assume that wave function obtained in this procedure 
is a good approximation for the true wave function, that 
is, it has a large overlap to true one. In practical shell 
model calculations, this assumption seems to be good 

El. 

Under this assumption, we consider how to construct 
a series of systematically approximated wave functions. 
The wave function of lowest minimum is written by 
|</5min}- The shell model orbits are grouped into upper 
and lower orbits, for instance, in the /p-shell, four or- 
bits are grouped by {/7/2} and {^3/2, /5/2,Pi/2}- By 
this grouping, we consider a following parameterized pro- 
jected wave functions: 

\^J,M (x)} =J29kPmk \^ix)) (3) 

K 

and 

\ipix))=l[aiix)\0) (4) 

where the creation operator ajj(x)'s are defined as 
'^q(^) — X] The corresponds to the 

i 

one of |(/3inin) and Xi takes 1 or x, depending on the 
group of i. The g's are determined at each x. When 
the X increases or decreases from unity, the angular mo- 
mentum projected energy of \ip (x)) increases because of 
\(p{x = 1)) = |</?min)- In this way, we can systematically 
generate a series of well- approximated wave functions 
near the lowest minimum continuously as a function of 



energy. This construction is not unique and we may find 
other methods. However this way to construct a series of 
wave functions works well for /p-shell calculations as we 
will show later. 

Note that, in general, we have no method to know 
whether or not a considered wave function is well- 
approximated to the true wave function unless informa- 
tion of the true state is available. Empirically an opti- 
mized projected wave function by variation after angu- 
lar momentum projection is known to describe deformed 
state well while we can not, strictly, know the precision 
of its approximation. In the present method, the energy 
variance can give an index of precision of approximation. 

For the extrapolation by these deformed bases, expec- 
tation value of (-ff^) is necessary. As we consider a real- 
istic description, we use a general two-body shell model 
effective interaction as, 

^ = ^eQclca+ ^ VapfScic'pCsCj (5) 
a a<p,f<6 

where e's and v's are single particle energies and two- 
body matrix elements, respectively. Then includes 
one-body, two-body, three-body and four-body terms. 
For instance, normal ordered form of four body term is 
shown by 

i?4 = ^ Vijklal3jSclcjclclcsCjCf3Ca (6) 

i<j<k<l,a<l3<'i<S 

where v can be given by summing up relevant w's and is 
totally antisymmetric for i,j,k,l and for a,(3,j,S. Ex- 
pectation value of this operator by angular momentum 
projected state can be carried out by Wick's theorem. 
We define generalized one-body density matrix as 

Pl3a = {^\ cicf3 IV')/ {v\ V') (7) 

where \ip) and lip) are deformed Slater determinants. For 
angular momentum projection, density matrices between 
different wave functions are necessary. By this one-body 
density matrix, we can define two-body density matrix 
as 

{ip\ clc^csc^ IV') 

PaP-iS = (^pV) ^ PiaPSf3 " P&aP-,fi- (») 

By the two-body density matrix, we can show the 
expectation value of the four-body term H/^ which can 
be given by six terms of product of two two-body 
density matrices. As we consider proton-neutron sys- 
tem and isospin projection is not introduced, two- 
protons and two-neutrons creation and annihilation part 
c]^ -c^ .,cJ, ^,cJ,_;Ci,_5Cj,_-),C7r,/3C^,Q, which is a dominant part 
of the (i?^) calculation, needs only one product of two 
two-body density matrices. Then such a normal ordered 
form of the 11^ operator is most efficient for numerical 
calculation. 

Next we consider how to reduce computational time 
by taking an advantage of time-reversal symmetry. By 
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solving eq. ^ ioi J = ground state, we can obtain the 
wave function with time-reversal symmetry. For such a 
case, it is known that an integral domain for angular 
momentum projection can be reduced. The projected 
matrix element of the Hamiltonian H^j^, is defined by 

^KK' ~ (V^l ^Prk' IV'): which is evaluated by a follow- 
ing integral: 

=w r^a [ dp r dw'K.K- m ■ (9) 

Jo Jo Jo 

The n stands for the Euler's angles. The ^, (Q) is a 
sum of relevant product of the Wigner's D-function and 
rotated matrix element (V'l HR{fl) \tp) where R{^1) is a 
rotating operator. Detailed form is written in Ref. [20| . 
for instance. Due to the reduced integral domain and 
relations of rotated matrix element concerning Euler's 
angles, computational time reduces 1/8 shorter. 

Furthermore, to reduce amount of computation, as a 
series of systematically approximated wave functions for 
the extrapolation of excited states, we use the same series 
of ground state. In this case, we can use the same rotated 
matrix elements for {H'^) calculation. Therefore we can 
simultaneously obtain extrapolations of several low-lying 
states in addition to the one of the ground state. 
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FIG. 1: Extrapolations of two kinds of extrapolation method 
are shown for the ground and first excited states of ®°Zn with 
the FPD6 interaction 16). The open circles are shown for 
energies and energy variances of t = 3 ~ 5 truncated spaces. 
Their extrapolations are shown by dashed lines. The filled 
squares and filled diamonds are energies and energy variances 
of deformed states with angular momentum projection for 
groupings (i) and (ii), respectively. The angular momentum 
projected energies and energy variances for J = and 2 from 
the Hartree Fock wave function are shown by open diamonds 
and are indicated by the word VBP. The energy of variation 
after angular momentum projection for J = state is labeled 
by the word VAP. 

In order to test the present extrapolation method, we 
take an example of ^°Zn in the /p shell model space, 



the dimension of which Hilbert space is largest in the 
/p shell model space and is 2,292,604,744 in M = 
space. We use the FPD6 effective interaction jl^. To 
test the present extrapolation, the exact ground state 
energy is needed. The exact ground state energy of this 
nucleus may be obtained by the state-of-the-art Lanczos 
shell model diagonalization if we take quite long compu- 
tational time. However here we use extrapolation method 
which we have proposed in Refs. The extrapola- 

tion method also use the scaling property of energy and 
energy variance while we combined it to Lanczos shell 
model diagonalization. By the minimum number t of the 
holes in the j-jj^ orbit, we can systematically define the 

truncation spaces as ^tf^^^^^r^^* (r = P3/2, /5/2,Pi/2)- 
Truncation spaces are labeled by t. As the t increases, 
its truncation space becomes larger and the lowest energy 
state |(/3f) of such a truncated space becomes a better ap- 
proximation to the true state. 

In Fig. 1, we use the truncation spaces with t =3, 4 
and 5. The dimensions of truncated shell model spaces 
are 4794006, 31077402 and 133533398 respectively. For 
each space, we diagonalize the shell model Hamiltonian 
matrix and obtain the approximated wave functions of 
the true ground state. Then we evaluate the matrix ele- 
ment of iP' . The extrapolation is shown in Fig. 1. For 
the diagonalization in truncated spaces, we use the shell 
model code MSHELLj^J. 

As deformed basis, Hartree-Fock energy is primarily 
important and is -243.03 MeV for s°Zn. First step of 
improvement of wave function is the angular momentum 
projection. The angular momentum projected energies 
for J = and 2 states from the Hartree Fock wave func- 
tion are shown in Fig. 1. Angular momentum projection 
improves energy and we can explicitly handle angular 
momentum quantum numbers. This approach is known 
to be the VBP (variation before projection). Next step 
is a variation after angular momentum projection. We 
can determine the minimum energy in the angular mo- 
mentum projected energy surface for J = 0. This energy 
is also shown in Fig. 1. By these energies and energy 
variances, we may extrapolate the true energy. In this 
case, this extrapolation is not so bad. However, rela- 
tion between these two wave functions is not clear and 
extrapolation from only two data is not reliable. 

For extrapolation, a series of systematically approxi- 
mated wave functions is important. To construct such 
a series, four /p shell orbits are grouped in following 
two ways: (i) {/7/2,P3/2} and {/5/2,Pi/2} and (ii) {/7/2} 
and {P3/2, /5/2,Pi/2}- For these groupings, we change x 
parameter in \Lp{x)). The dependence of x in ajj(a;) is 
simple, but for the angular momentum projected energy, 
the dependence of x is no more simple and is non-linear. 
However by rescaling the energy as a function of energy 
variance, we can obtain a simple scaling again. In Fig. 
1, we take a;=1.0, 0.8, 0.75, 0.7, 0.65 and 0.6 for group- 
ing (i) and a;=1.0, 0.8, 0.7, 0.6, and 0.5 for grouping (ii). 
These data of energy to the energy variance for ground 
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and excited states are well-aligned and we can find fine 
proportionalities, by which we can make sure extrapola- 
tion. The energies of the present extrapolation method 
agree with those of the previous method within 0.1 MeV. 
In the respect of the computational amount, previous ex- 
trapolation with spherical basis [T^ needs four days 
while this new approach needs less than half day in this 
calculation. Therefore the present extrapolation method 
is superior to the previous one. 

In the present method, we solve eq. Q for J = 
ground state while extrapolations are carried out for sev- 
eral low-lying states from the series of wave functions 
for J = state. This method has an advantage for 
drastically reducing computational time. However for 
higher spin states, the approximation may become worse. 
Therefore, we test what extent the extrapolation works 
well for several low-lying yrast states by this approach. 
Here we take ^®Cr as an example. For comparison, we 
can give exact energies for yrast states by Lanczos di- 
agonalization. In the left panel of Fig. 2, we show the 
exact level scheme. In the right panel, the energies as a 
function of energy variances for x=1.0, 1.1, 1.2 and 1.3 
are plotted. Four fp shell orbits are grouped in a follow- 
ing way: 1/7/2} and {P3/2, f5/2,Pi/2}- Because x = 1.0 
corresponds to the minimum energy, energies with a; > 1 
are higher than minimum energy. By this construction, 
for a series of wave functions x > 1 (a: < 1), excitation 
from the f'j/2 orbit to other upper three orbits increases 
(decreases). As more excitation is needed, in general, to 
generate high spin states, we consider a series for x > 1. 
In Fig. 2, the extrapolated energies agree with the ex- 
act ones within 0.1 MeV, which means that a series of 
wave functions for J = state works for other low-lying 
states. However, minimum of variance for higher spin 
state increases. This means deterioration of approxima- 
tion, while such deterioration is somewhat compensated 
by the extrapolation procedure. Thus the present ap- 
proach works well for low spin yrast states. 

For higher spin, deterioration of approximation is evi- 
dent. This point is expected to be improved by optimiz- 
ing a deformed wave function for each spin, which needs 
longer computational time. Moreover, here we show en- 
ergies for J < 8 states while this nucleus is known to have 
a backbending [22| . For such a large change of structure 
in wave functions as a function of spin, such an optimiza- 
tion becomes necessary. 

Finally we discuss an advantage of the present ex- 
trapolation method to other extrapolation methods. In 
an open-shell nuclei, deformed states often appear due 
to residual interaction. To describe deformed states by 
spherical bases, we need almost entire spherical bases. 
Consequently truncation is not so effective and the con- 
vergence of energy as a function of the number of spher- 
ical basis state becomes very slow. According to the 
ECM0, necessary number of spherical basis state in- 
creases exponentially as we need more precise energy. 
On the other hand, it is evident that deformed basis can 
handle deformed solution more naturally. In the present 



-25 



-30 



- 


^C) 


■ g+ 








" 6* 




- 4*_ 




- 2*_ 
. 0"- 


^^^^ ""Cr (KB3 interaction) 



0.02 0.04 

Energy variance 



FIG. 2: Left panel: The exact level scheme is shown. Right 
panel: Extrapolations of the present method are shown for 
several low- lying excited states (0"*", 2^, 4"*", 6^ and 8^ ) of 
**Cr with the KB3 interaction |23l . The angular momentum 
projected energies and energy variances from the Hartree Fock 
wave function are shown by open diamonds. 



paper, we use always one Slater determinant. This fact is 
in sharp contrast with other extrapolation methods and 
becomes a distinct advantage for larger shell model prob- 
lems, where a physically relevant truncation by spher- 
ical basis becomes severely difficult. Another merit of 
the present method is to handle all nucleus in the con- 
sidered model space within almost the same amount of 
computation. In extrapolation methods based on the 
spherical shell model, as the number of valence nucle- 
ons reaches half of degeneracy of orbits, amount of com- 
putation rapidly increases. Therefore feasibility of the 
present method is wider than other extrapolation meth- 
ods. Moreover the present method has a close link to 
mean-field or projected mean-field theory. It is also an 
advantage that physical interpretations of the obtained 
results based on mean-field is easy. 

In summary, we have proposed an extrapolation 
method in shell model calculations by deformed wave 
functions with angular momentum projection. Our ex- 
trapolation technique relies on a scaling property be- 
tween energy and energy variance of a series of systemat- 
ically approximated wave functions. To apply it to shell 
model calculations, a key issue is how to construct such 
a series. In our previous papers |l2l H^ . we consider 
physically relevant truncation subspaces and we solve 
shell model problem within these subspaces by Lanczos 
method. Then we obtain such a series and apply the ex- 
trapolation technique. On the other hand, in the present 
paper, we consider a deformed wave functions with angu- 
lar momentum projection. By variation-after-projection, 
we can solve shell model problem under a restriction of 
single Slater determinant with angular momentum pro- 
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jection. This wave function often becomes a relevant ap- 
proximation to the true state, especially when nucleus is 
deformed. We have shown that a series of approximated 
wave functions can be generated by changing their struc- 
ture systematically. To evaluate the energy variances 
by such an angular momentum projected Slater deter- 
minant, the Wick's theorem is used. As we do not use 
isospin projection for proton-neutron systems, most time 
consuming part of {H"^) calculations is rather simplified. 
Time reversal symmetry is also useful to reduce the nu- 
merical calculations concerning angular momentum pro- 
jection. By these developments, we have shown that the 
extrapolation method with deformed basis works well. 

In the present paper, for minimizing necessary amount 
of computation, we introduced several restrictions. How- 
ever, for more precise extrapolation and description of 
excited states with the same quantum numbers, we can 
introduce a diagonalizaiton with several optimized Slater 
determinants and optimization of deformed bases for 
each spin. Such important extensions will be able to be 
realized with a help of parallel computation. 



As applications of the present method, fpg shell is 
important, whose orbits are /5/21 P3/2 and gg/2- 

The ground state of A ~ 80 nuclei is determined by 
a subtle competition between oblate and prolate de- 
formed states and occupation of the /2 orbit is essen- 
tial, which means that relevant truncation spaces become 
quite huge. On the other hand, deformed states can be 
well-approximated by well-optimized Slater determinant 
with angular momentum projection. The present extrap- 
olation method can handle such a subtle competition. 
Moreover, the present extrapolation may also apply to 
the study of fractional quantum Hall effect formulated 
on the Haldane sphere 0|. Such studies are being in 
progress. 
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